On the existence, uniqueness and nature of Caratheodory and Filippov 
solutions for bimodal piecewise affine dynamical systems 



L. Q. Thuan a - b , M. K. CamlibeP 



a Johann Bernoulli Institute for Mathematics and Computer Science, University of Groningen, P.O. Box 800, 9700 AV 

Groningen, The Netherlands 
b Department of Mathematics, Quy Nhon University, 170 An Duong Vuong, Quy Nhon, Binh Dinh, Vietnam 
c Department of Electronics and Communication Engineering, Dogus University, Kadikoy 34722, Istanbul, Turkey 



Abstract 

In this paper, we deal with the well-posedness (in the sense of existence and uniqueness of solutions) 
and nature of solutions for discontinuous bimodal piecewise affine systems in a differential inclusion setting. 
First, we show that the conditions guaranteeing uniqueness of Filippov solutions in the context of general 
differential inclusions are quite restrictive when applied to bimodal piecewise affine systems. Later, we 
present a set of necessary and sufficient conditions for uniqueness of Filippov solutions for bimodal piecewise 
affine systems. We also study the so-called Zeno behavior (possibility of infinitely many switchings within 
a finite time interval) for Filippov solutions. 

Key words: Piecewise affine systems, well-posedness, existence and uniqueness of solutions, Caratheodory 
solutions, Filippov solutions, one-sided Lipschitz condition. 



1. Introduction 

A piecewise affine dynamical system is a special type of finite-dimensional, nonlinear input-state-output 
systems with the distinguishing feature that the functions representing the system of differential equations 
and output equations are piecewise affine functions. Such systems arise in various contexts of system and 
control theory such as variable structure systems [15] . bang-bang control [12], and linear relay systems [HI2U]. 

Piecewise affine functions which describe the dynamics of a piecewise affine dynamical system are not 
necessarily continuous. As such, piecewise affine dynamical systems form a subclass of discontinuous dynam- 
ical systems (see [5] as an excellent survey). An immediate consequence of discontinuous dynamics is that 
the existing results of mainstream smooth nonlinear systems and control theory (see e.g. [16) ) cannot be 
indiscriminately applied to piecewise affine dynamical systems. This departure from smooth systems begins 
already from the definition of a notion of solution. Indeed, meaning of a solution of a differential equation 
given by continuous functions is rather straightforward whereas it becomes a much more complicated matter 
in the absence of continuity. 

The typical framework to deal with discontinuous dynamical systems is the framework of differential 
inclusions (see e.g. |18j). Roughly speaking, one replaces a differential equation with discontinuous right- 
hand side (see e.g. Filippov's seminal work [2]) by a differential inclusion given by a set- valued mapping. 
There are several ways of defining a set-valued mapping (and hence a differential inclusion) for a differential 
equation with discontinuous right-hand side. Each of these leads to a different solution concept (see e.g. [HE] 
for details) such as Caratheodory, Krasovskii, Filippov, and Euler solutions. A good deal of the literature 
on differential inclusions is devoted to the investigation of existence and uniqueness of solutions. Typically, 
existence of solutions is guaranteed by less restrictive conditions than those for uniqueness. 
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In this paper, we focus on a particular class of piecewise affine dynamical systems, namely bimodal 
piecewise affine systems without external inputs. The main goal of the paper is to investigate existence, 
uniqueness, and nature of solutions (in the sense Caratheodory and Filippov) for this class of systems. 
It turns out that existence of Filippov solutions immediately follows from the existing results for general 
differential inclusions. However, existing conditions ensuring uniqueness for general differential inclusions 



are quite restrictive in the context of piecewise affine dynamical systems (see Theorem 2.6 in Section Yl 



Motivated by this fact, we turn our attention to tailor-made conditions for bimodal systems. The main 
results of this paper are a set of necessary and a set of sufficient conditions for uniqueness of Filippov 
solutions of bimodal systems. Furthermore, these results provide necessary and sufficient conditions for the 
existence and uniqueness of Filippov solutions for bimodal piecewise linear systems. 

A curious phenomenon in the context of discontinuous dynamical systems is the so-called Zeno behavior 
(see e.g. (101 lllj ) which refers to infinitely many switchings in a finite time interval. Presence of such 
behavior causes serious difficulties not only in analysis and design but also in simulation of such systems. As 
a by-product of our main results, we obtain conditions under which Zeno behavior is ruled out for bimodal 
piecewise affine systems. 

Well-posedness of piecewise affine dynamical systems has received considerable attention in the last 
two decades. In jS], the authors consider bimodal piecewise linear systems. They work with what we call 
forward Caratheodory solutions and provide necessary and sufficient conditions for existence and uniqueness 
of these solutions. Forward Caratheodory solutions rule out the possibility of left accumulation points for 
switching instance by their very definition. In this paper, we consider more general bimodal systems, namely 
bimodal piecewise affine systems. Also we work not only with forward Caratheodory solutions but also with 
Filippov solutions. As such, the main result of [5] becomes a special case of our main results. In 20], well- 
posedness of linear relay systems was addressed for forward Caratheodory solutions and sufficient conditions 
for uniqueness were presented. A linear relay system with a single relay boils down to a bimodal piecewise 
affine system as studied in this paper (see Example 2.1 ). The results presented in this paper show that the 
very same conditions of |20j ensure uniqueness of Filippov solutions for this case. The paper [3] studied 
linear relay systems with a single relay and provided sufficient conditions for the uniqueness of Filippov 
solutions. Also the results of @] can be recovered from our main results. Another related paper is [5] which 
considers Filippov solutions bimodal piecewise linear systems. The results of [5] can also be recovered as a 
special case from our main results (see Corolloary 3.5 1. 

Zeno behavior of systems that are closely related to piecewise affine dynamical systems has been consid- 
ered in [9 for a class of linear complementarity systems, in [15] for conewise linear systems, in |17j for linear 
relay systems with a single relay, and in [14] for continuous bimodal piecewise affine systems. 

The paper is organized as follows. In Section [2] we introduce the object of the study in this paper, i.e. 
bimodal piecewise affine dynamical systems in a differential inclusion setting. We also present two examples 
of such systems and define what a solution means for such systems. This is followed by a discussion on 
existence of solutions as well as a discussion of the restrictiveness of conditions that guarantee uniqueness 
of Filippov solution of general differential inclusions when applied to bimodal systems. Section [3] presents 
the main results related to the uniqueness of Filippov solutions whereas we investigate Zeno behavior of 
bimodal systems in Section [4] In Section [5j we present the proofs of the main results. Finally, the paper 
closes with the conclusions in Section |6l 



2. Bimodal piecewise affine systems 

Consider the differential inclusion 



with 



x(t) e F(x(t)) (1) 



'{Aix + eJ if c T x + f < 

F(x) = { {Aix + e 1 ,A 2 x + e 2 } if c T x + / = 
{A 2 x + e 2 } if c T x + f > 
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where x € W 1 , A%, A 2 € M. nxn , e\, 62, c € IR n and /gt. Also consider the convexified differential inclusion 

x(t) e G(z(t)) (2) 

with 

{{Aix + ei} if c T x + / < 

conv({Aix + ei, A 2 x + e 2 }) if c T x + f = 
{A 2 x + e 2 } if c T .x + / > 

where conv stands for the convex hull. 

Throughout the paper, we call the systems of the form ([!]) and |2| bimodal piecewise affine systems. In 
the sequel, we investigate existence and uniqueness of different kinds of solutions of bimodal piecewise affine 
systems. 

Before elaborating on the solution concepts for these systems, we provide some examples of bimodal 
piecewise affine systems. 

The first class of examples consists of linear systems with ideal relay elements which serve as an idealized 
models of Coulomb friction, bang-bang control, etc. 

Example 2.1. (Linear relay systems [HHU]) Consider the linear relay system 

x(t) = Ax(t) + bu{t) 
y(t) = c T x{t) 
u{t) e - B gn(y(i)) 

where sgn is the set-valued relay function defined by 

({-1} if y < 
sgn(y) = l[-l, 1] if y = 
[{1} if y > 0. 

Clearly, such a linear relay system is a bimodal piecewise affine system of the following form 

( {Ax(t) + b} if c T x(t) < 

±(t) e < conv{{Ax{t) + b, Ax(t) - b}) if c T x(t) = 
(_{Ar(t)-6} if c T x{t) >0. 

Note that linear relay systems are particular cases of ^ where A\ — A2 and / = 0. Well-posedness 
of such systems were studied in [4] and [20] , Another example of bimodal piecewise affine systems is the 
following water tank system. 

Example 2.2. Consider the two tank system depicted in Figure [T] 

The deviations of the water level from the bottom of the first and second tank are denoted by x\ and 
x 2 , respectively. Let u be the constant flow of water into the first tank. The valve V is opened if a; 2 < 1 and 
closed if xi ^ 1. By defining the state x = col(xi, x-i) and taking all involved parameters unity, one obtains 
the following equations describing the dynamics of the system 

if x 2 - 1 < 





-1 


" 




u 


x = 


1 


-1 


x + 








-1 " 




"0" 


X = 


1 -1 


X + 






]£x 2 -l>0 



As such, this system is of the form of a bimodal piecewise affine system (|T 
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Figure 1: Water level regulator 



We now turn our attention to formalizing what will be meant by a solution of the system |T]). There 
are many ways of defining a solution for a differential inclusion (see e.g. [HE]). In this paper, we focus on 
Caratheodory and Filippov solution concepts. 

Definition 2.3. An absolutely continuous function x : K — > K" is said to be a solution of the system ([!]) 
for the initial state £ in the sense of 

• Caratheodory if x(0) — £ and x satisfies the differential inclusion ([!]) for almost all (el. 

• forward Caratheodory if it is a solution in the sense of Caratheodory, and for each if there exist 
i € {1,2} and e t * > such that 

x{t) = A lX {t) + e< and [c T x{t] + f] < (3) 

for all te(t*,t*+ e t »). 

• backward Caratheodory if it is a solution in the sense of Caratheodory, and for each t* there exist 
i € {1,2} and et* > such that 

x{t) = Aix(t) + ei and [c T x{t) + /] < (4) 

for all t e (t* - e t .,t*). 

• Filippov if ir(0) = £ and x satisfies the convexified differential inclusion ^ for almost all t G R. 

Clearly, every Caratheodory solution is a Filippov solution since F{x) C G(x) for all x £ W\ However, 
not every Filippov solution is a Caratheodory solution in general. 

When the right hand side of (JlJ is single- valued and hence is Lipschitz continuous, that is the implication 

c T x + f = =*> A x x + ei = A 2 x + e 2 (5) 

holds, Caratheodory and Filippov solutions coincide. In this case, existence and uniqueness of solutions are 
guaranteed by the theory of ordinary differential equations. 

In general, existence of solutions of the differential inclusion ([TJ readily follows from the theory of 
differential inclusions (see e.g. [2j Theorem 2.7.1]). 

Proposition 2.4. There exists a solution of the differential inclusion |TJ in the sense of Filippov for each 
initial state. 
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In the sequel we focus on the uniqueness of Filippov solutions for the system . 

Definition 2.5. We say that a Filippov solution for the initial state £ is right-unique (left-unique) if for 
any Filippov solution x' for the initial state £ there exists e > such that x(t) = x'(t) for all t G [0,e) 

(te(-e,o]). 

The main goal of the paper is to present necessary and/or sufficient conditions for uniqueness of Filippov 
solutions that are tailored to bimodal piecewise afhne systems of the form ([!]). To motivate these new 
conditions, we hrst review one of the most typical uniqueness conditions that is employed in the literature 
of (general) differential inclusions and discuss its limitations for bimodal piecewise affine systems. 

A set- valued mapping H : R" =4 W l is said to be one-sided Lipschitz (see e.g. [B]) if there exists a 
number L such that 

{xi - x 2 ) T {yi - y 2 ) «S L\\xi ~ X2W 2 (6) 

for all x\, x 2 belonging to the domain of H, y x G H(x\) and y 2 G H.{x%). 

The following theorem presents necessary and sufficient conditions for one-sided Lipschitzness of the 
set- valued mapping G. 

Theorem 2.6. The set-valued mapping G is one-sided Lipschitz if and only if there exist a vector g G R™ 
and a number pi ^ such that 

Ai — A 2 — gc T and e x — e 2 = fg + [ic. 

Proof. For the 'if part, suppose that we have A\ — A 2 = gc T , e\ — e 2 = fg + /ic where /i ^ 0. Then, the 
piecewise affine function which is defined by 



G(x) 



Aix + ei - y if c T x + f sc 
A 2 x + e 2 + y if c T x + f ^ 



is continuous, so that it is globally Lipschitz continuous (see [7 ). Observe that G(x) = G(x) + [1 — 2A(x)] — 
where 

( {0} if c T x + f < 
X(x) = I [0, 1] if c T x + f = 
[{1} if c T x + f> 0. 

Thus, for any Xi G K. n and j/j G G(xi), there exists y~i G G(xi) such that yi = yi + (1 — 2Xi)^c where 
Ai G X(xi), and then 

(xi - x 2 ) T (yi - y 2 ) = (xi - x 2 ) T [{y x - y 2 ) - (j,(X 1 - A 2 )c] . 

This together with the observation that // ^ and (Ai — X 2 )(xi — x 2 ) T c implies that 

(xi - x 2 ) T (y! - y 2 ) < (xi - x 2 ) T {yi - y 2 ). 

On the other hand, the Cauchy-Schwarz inequality and Lipschitzness of G implies (x\ — x 2 ) T (yi — y 2 ) ^ 
L\\xi — x 2 \\ 2 where L denotes a Lipschitz constant of G. Thus, we have 

(xi - x 2 ) T (y 1 - y 2 ) < L\\xi - x 2 \\ 2 

for all Xi,x 2 G R" and t/j G G(xi), and hence G is one-sided Lipschitz. 

For the 'only if part, we suppose that G is one-sided Lipschitz. Let £_ = {x | c T x + f ^ 0} and 
S + = {a; I c T x + f ^ 0}. Let G E_, x 2 G E+, and let x be such that c T 5 + / = 0. For a G (0, 1], define 
x\ = ax \ + (1 — a)x and x' 2 = ax 2 + (1 — a)x. Clearly, x[ G E_ and x' 2 G E + . Since G is one-sided Lipschitz, 
one has 

\{A x x' x + e x ) - (A 2 4 + e 2 )f{x' 1 - x' 2 ) ^ L\\x' x - x' 2 \\ 2 , 



or equivalently 
-(1-a) 



Q 



{(A 1 - A 2 )x + (e x - e 2 )} + {(Am + ei) - (A 2 x 2 

A 2 )x + ( ei - 



e 2 )}] (xi - 2:2) < L\\xi - x 2 \ 



e 2 )] T (a;i — x 2 ) for all xi € S_ and 

(7) 



By taking sufficiently small a, we obtain that [(A\ 
X2 G £+. This implies that 

(Ai - A 2 )sc + (ei - e a ) G (£_ 

for any S with c T x + / = 0. Here the notation S° denotes the polar cone of the set S that is S° — {y \ 
x T y 0, Vie S}. Then, we get 

{A x - A 2 )(kerc T ) + (Aj - A 2 )x + (e x - e 2 ) C (S_ - £+)° = {ac \ a ^ 0} 

for fixed a; satisfying c T x + f = 0. Since the left hand side is an affine set and the right hand side is a cone, we 
can conclude that {A\ — A 2 )(kerc T ) = {0}. Hence, A\ — A 2 = gc T for some g. Then it follows from |7]) that 



ei-e 2 -fge 
This means that e\ — e 2 



Note that £_ 



fg + fie for some /1 ^ 0. 



{x I c 2 x < 0}. Hence, we have (£_ — £+)° 



{ac I a ^ 0}. 



Theorem |2.6| shows that we can employ one-sided Lipschitzness in order to conclude uniqueness of 
Filippov solutions of bimodal systems only under quite restrictive conditions. Note that these conditions are 
met for a linear relay system as in Example |2 . 1 1 only if b = ac for a nonnegative real number a. Note also 



that the conditions of Theorem 2.6 are never met for a two-tank system as in Example 2.2 Motivated by 
the restrictiveness of one-sided Lipschitzness, we investigate tailor-made uniqueness conditions for bimodal 
systems in the sequel. 

3. Uniqueness of solutions 

In this section, a set of necessary and a set of sufficient conditions for right-uniqueness of solutions 
of the bimodal piecewise affine systems ([l]) will be provided. These conditions are less restrictive than the 
conditions guaranteeing the one-sided Lipschitz condition. To do so, we need to introduce some nomenclature 
as follows. For a vector v, we write v >- if it is nonzero and the first nonzero entry is positive. We write 
v >z if either v = or v >~ 0. Similarly, we write v -< when — v >~ and v ■< when — v >z 0. The 
observability index of the pair (c T ,Ai) is denoted by hi, that is the largest integer such that the matrix 



col(c 
such that 



A, 



,c T Ah 



has full row rank. Note that for each k ^ 1 there exists a unique matrix Pf £ 



tkxhi 



\c T A'l' +1 - 




r ° T 1 


c T A h>+2 


1 i 


c T Ai 




_ c T A h i+ k_ 




c T At\ 



(8) 



We now present the main results concerning the well-posedness of solutions of the bimodal piecewise 
affine system ([TJ in the three theorems below. 

Theorem 3.1. Let h := min{hx, h^}- Consider the statements: 

1. Every Filippov solution of ([T]) is right-unique. 

2. Every Filippov solution of ([T]) is both a forward and backward Caratheodory solution. 

3. There exist an integer k with 1 k h + 1 and a (k + 1) x (fc + 1) lower triangular matrix M with 
positive diagonal elements such that 



- C T - 

c T A 1 


= M 


r c T ' 

c T A 2 


5 


" / " 

T 

c ei 


y M 


r / " 

T 

c e 2 


c T A\_ 




?A\_ 








_c T A\^e 2 _ 



G 



4- There exists a (h + 1) x (/i + 1) lower triangular matrix M with positive diagonal elements such that 



- C T - 

c T A 1 


= M 


\ c T - 

c T A 2 


1 


\ f ' 

T 

c e\ 


= M 


\ f 1 

T 

c e 2 


c T A\_ 




fA\_ 




fA\~ x e x _ 




c T A\- x e 2 _ 



and either h\ < h 2 and (FA^ei — p{ e" 1 > or hi > h 2 and c 1 A 1 ^ 2 e 2 — p 2 e^ 2 < where pi is uniquely 
determined from 

c T A^ +1 =pfT l h \i = 1,2. 

5. The observability indices h\ and h 2 are the same and there exists a (h + 2) x (h + 2) lower triangular 
matrix M with positive diagonal elements such that 



T Aha 



- C T - 

c T A 1 


= M 


\ c T - 

c T A 2 


1 


1 

T 

c e\ 


= M 


f 

T 

c e 2 


c T A h +\ 




c T A h +\ 




c T A h iei _ 




c T A h 2 e 2 _ 



6. The following implication holds 



r c T - 




/ 


c T A 1 




T 

c ei 








c T A h + l _ 




c T A\e x _ 



= O^AiC + ei = A 2 e, + e 2 . 



Then, the following implications hold: 

A. [l]^ (d org or[5} 

B. and |5| =^H 

C. (|]and[6) 

D. (|5]and[6| ^[T] 

A proof of this theorem will be presented in Section [5| Note that this theorem provides only a set of 
necessary and a set of sufficient conditions, but not necessary and sufficient conditions in general. The 
following example (see [H Eq's. (13) and (14)]) illustrates the gap between the necessary and the sufficient 
conditions. 

Example 3.2. Consider the bimodal piecewise affine system ([I]) with c T = [l 0] , / = 0, 



Ai = A 2 



1 
1 




, and ei 



-e 2 



For this system, it can be verified that the third statement of Theorem 3.1 is satisfied with k = 3. However, 
as it has been shown in [4], there are infinitely many Filippov solutions for the zero initial state. 



However, the third statement of Theorem |3.1| with k = 1 is sufficient for right-uniqueness of Filippov 
solutions as stated in the following. 
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Theorem 3.3. If there exists a 2 x 2 lower triangular matrix M with positive diagonal elements such that 



c 

c T A 1 



= M 



T 

c e\ 



y M 



T 

c 1 e 2 



(9) 



then right-uniqueness of Filippov solutions holds at any state of the state space K. n . 

Proof. Due to (foh , for each state £ satisfying c T £ + / = at least one of the inequalities c T 'Ai£ + c T e\ > 
or c T A 2 £ + c T e 2 < is satisfied. By [21 Theorem 2.10.2], every Filippov solution is right-unique. ■ 



" c T ' 




' c T ' 




/ 




f 


c T A 1 


= M 


c T A 2 




T 

c e\ 


y M 


T 

c 1 e 2 


c T A\ 




c T A\ 




c T A 1 e 1 




_c T A 2 e 2 _ 



Further, the third statement of Theorem |3.1| with k = 2 is sufficient for right-uniqueness of Filippov 
solutions for some initial states. 

Theorem 3.4. // there exists 3x3 lower triangular matrix M with positive diagonal elements such that 

(10) 

then for the system ([T]) 7 right-uniqueness of Filippov solutions holds at all states from K"\f2 where Q is the 
set of all £ € R™ such that 

r cT 1 I" / 1 

Ij e + c T e 3 = and {-l) k {c T A 2 k i + c T A k e k ) < 
tf\ [_ cT Ajej_ 

for some k, j G {1, 2}, k ^ j. 

A proof of this theorem will be given in Section [5j 



Theorem |3.3| and Theorem |3.4| present two particular cases under which the third statement of The- 
orem |3.1| becomes sufficient as well as necessary for right-uniqueness of Filippov solutions for bimodal 
systems. Another interesting particular case occurs when there are no affine terms in the dynamics, that is 
when e\ = e 2 — and / = 0. In this case, one can state necessary and sufficient conditions (see also [5]) as 
in the following. 

Corollary 3.5. Consider the system |l]) with e\ = e 2 = and / = 0. Then, every Filippov solution of ([I]) 
is unique if and only if the following statements hold: 



1. hi = h 2 . 

2. There exists an (hi + 1) x (hi 
that 



1) lower triangular matrix M with positive diagonal elements such 



r c T - 




r ° T i 


c T A 1 




c T A 2 




= M 




c T A\\ 




fA\\ 



3. The following implication holds 



c T 
c T Ai 



c T A\ 1 



Aix = A 2 x. 



Proof. Note that the third and fourth statements of Theorem 3.1 never holds as ei = e<z = and / = 0. 
Then, the first statement of Theorem |3 . 1| holds if and only if the fifth and the sixth hold. ■ 
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4. Switching behavior 



In this section, we investigate mode switching behavior of bimodal systems. We say that a time instant 
t* G K is a non- switching time for a Filippov solution x if there exist an interval (t* — e, t* + e) and an index 
i with i € {1,2} such that 



x(t) = Aix(t) + e { and (-1)' 



x(t) + f] < 



for all £ G (t* — e, i* + e). We say that a time instant t* G K is a switching time for a Filippov solution a; if 
t* is not a non-switching time for the same solution. 

The distribution of the switching times along a solution is an important issue for various reasons. For 
instance, the so-called event-driven simulation methods (see e.g. [3]) may fail if the switching times accu- 
mulate around a point. This type of phenomenon is known as Zeno behavior in hybrid systems literature. 
We say that a time instant t* G R is a left/right Zeno time for a Filippov solution x if for each e > the 
interval (t* ,t* + e)/(t* — e,t*) contains a switching time for the same solution. 

A Fillipov solution will be called left (right) Zeno-free if there exists no left (right) Zeno time for it. 
Also, we say that the system ([I]) is left (right) Zeno-free if all its Fillipov solutions are left (right) Zeno-free. 
Further, we say that the system |l]) is Zeno-free if it is both is left and right Zeno-free. 

With these preparations, we can state the following sufficient condition for Zeno-freeness. 

Theorem 4.1. Suppose that the observability indices h\ and hi are the same. Let h = hi = hi. If there 
exists a (h + 2) x (h + 2) lower triangular matrix M with positive diagonal elements such that 



r c T - 

c T A 1 


= M 


r c T - 

c T A 2 




/ 

T 

c e\ 


= M 


f 

T 

c e 2 


c T A h 1 +\ 




c T A h 2 +\ 




_c T A\e x _ 




c T A>ie 2 _ 



and the implication 



I c T - 




f 1 


c T A 1 




T 

c ei 








c T A h 1 +\ 




c T A h iei _ 



= 0^Ai£ + ei =A 2 £,+e 2 



holds, then the system ([T]) is Zeno-free. 

Proof. Note that the system ([lj is Zeno-free if and only if any Filippov solution is both forward and 
backward Caratheodory solution. As such, the claim follows from Theorem |3.1|C| ■ 



5. Proofs 



5. 1 . Proof of Theorem 3. 1 



To prove Theorem |3.1| we first consider affine systems and introduce some notations. An affine system 
E(A,e,c T J) is given by 



x = Ax + e 
y = c T x + f 



(11a) 
(lib) 



where x G K™ is the state, y G K is the output, and all involved matrices are of appropriate sizes. By x(t; £) 
and y(t;£), we denote the state and the output of the system (11) for the initial state £, respectively. We 
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define the sets 



We 



= {£ | 3e > such that y(t; £) < 0, Vt G (0, e)}, 
= {£ | 3e > such that y(t; £) = 0, W G (0, e)}, 
= {£ | 3e > such that y(t;£) > 0, V£ G (0,e)}. 



Since y(t;^) is a real-analytic function for each initial state £, its sign on a small time interval (0, e) is 
completely determined by the values of its derivatives at t = 0, i.e. j/ fc )(0;£) for ft = 0, 1, .... Note that 
y(°\0;£) = c T £ + / and i/ fc )(0;£) = c T ^4 fc £ + c T A k ~ 1 e for ft > 1. Let i/ denote the observability index 
of the pair (c T , A). It can be verified that c T A 1 ^ + c T A l ~ 1 e = for all I with v + 1 ^ £ ^ 1 implies 
c T A f ^ + c T A e = for all I ^ 1. Together with the analyticity of the output y(t; £), this observation leads 
to the following immediate characterizations of the W-sets: 




{< 



- C T - 




/ 


c T A 




T 

c e 








c T A»+\ 




_c T A"e_ 




(12) 



Note also that 



Wvuw?u w^ 



u ffj U fV^ — jR". (13) 

In the sequel, we often use the W-sets corresponding to the two modes of the system ([lj . For brevity, we 
denote = Sj(^4j, e,, c T , /) and define 



W 2 °:=W£ 2 W+:=W+ W 2 := Wg 2 U W+ . 



(14a) 
(14b) 



We also define 



rpk 



c T Ai 



s T A fc 



and ey 



T 

c e< 



c A 



for ft ^ 1. With these preparations, we are ready to prove Theorem 3.1 



5.1.1. Proof of Theorem \3J\Z\ 

Right-uniqueness of every Filippov solution of the system ([I]) implies that 

w x ~ n w 2 = and Wi n w^ = 0. 



In view of (12 1 and (14), therefore, we have 

Jn+i 



J»a+1 



-< o 

> 



y o 



(15) 

(16a) 
(16b) 



as necessary conditions for right-uniqueness of every Filippov solution. Note that a consequence of the first 
implication is that 

Ti^ + e^O T^ + e§=<0. (17) 

In order to formulate this condition in terms of the parameters of the system ([I]), we invoke the following 
lemma which was proven in [13 . 

Lemma 5.1. Let Pi be an m x n matrix of full row rank and qi be an m-vector for i = 1,2. Then, the 
following statements are equivalent: 
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1. P\X -< q\ implies P 2 x -< q 2 . 

2. P\x < qi implies P2X -< q 2 . 

3. Either 



Pi = MP 2 and qi = Mq 2 



for some m x m lower triangular matrix M with positive diagonal elements, or there exist £ with 
1 £ ^ m and £ x ^ lower triangular matrix M with positive diagonal elements such that 

pf ] =A/pM and qf<Mq [ ? 

where the notation mV°' denotes the first k rows of a matrix/vector. 



Note that T/ 1 is of full row rank for i = 1,2 since h = nrin{/ii, h 2 }- Then, it follows from (17) and 
Lemma 15. II that either 



T h 



MT' 2 l and e 



Me 



(18) 



for some (h + 1) x (h + 1) lower triangular matrix M with positive diagonal elements, or there exist £ with 
1 ^ £ ^ h + 1 and £ x £ lower triangular matrix M with positive diagonal elements such that 



T, 



and ef _1 >- Afe^ 1 . 



(19) 



Suppose first that the latter holds. Since T-j 3 = c T = T° and ej = / = e^, the case £ = 1 is impossible. 
Then, we get 

T* = MT 2 fc and ej >- Me\ 

for some A; with 1 k /i and (fc + 1) X (k + 1) lower triangular matrix M with positive diagonal elements. 
This is nothing but the statement [3j Thus, it remains to show that ( 18 1 also implies the statement [3] or [4] 
or [5] To do so, we consider two cases: hi = h 2 and hi ^ h 2 . 

For the case h\ ^ h 2 , we prove that (18) implies statement [4] Note first that the first part of statement 
[1] holds due to (18). For the second part, we only need to prove for hi < h 2 . For the case hi > h 2 , the proof 
is similar. 

Since hi < h 2 and T% 2 is of full row rank, there exists £ such that T^f + eJj 1 = and T% 2 £ + e£ 2 y 0. 
Together with (16b) and (18), this implies 

Ii hl £ + ef 1 = and Z? 1+1 f + ei H+1 >- 0. 

This immediately implies c T A^ei — p^e^ 1 > 0. 

For the case hi = h 2 = h, we prove that (18) implies the statement [3] or [5] To this end, first note that 
it suffices to show that either 



rph + l 



MT% +1 



or 



T^ +1 = MTz +1 



and 



and 



h+l 



y Me. 



h+l 



=1 = Me 2 



h+l 



(20) 



(21) 



for some (h + 2) x (h + 2) lower triangular matrix M with positive diagonal elements. Indeed, the former 
would imply the statement [3] and the latter [5] In what follows, we will construct a matrix M which will 
satisfy one of these two conditions. To do so, note that 



c^+^pfTf with i= 1,2 
for some Pi,p 2 £ R ft+1 as h is the observability index of both (c r , and (c T , A2). Define 



(22) 



M 

T 

q a 
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h+1 



MT 2 +1 . Therefore, it remains to show that we can 



for some a. It follows from (TL8b and fl22) that T{ 
choose a > such that either e^ +1 >- Me' 2 l+1 or e'^ 1 = Afe£ +1 . Since ef = Afe^, e^ +1 >- Afe 2 l+1 holds if 
and only if c T Afei > g T e§ + acF A\e 2 and e^ +1 = Me 2 +1 holds if and only if c T A\e x = q T e' 2 l + ac T A 2 l e 2 . 
As such, it is enough to show that we can choose a > such that A\e\ > q 7 ' e 2 + ac T A 2 e 2 . By using the 
definition of q, we see that the last inequality is equivalent to 



(c T A h iei -pK) > a{c T A h 2 e 2 -p T 2 e 



T 4) 



(23) 



Since Tf is of full column rank, there exists £o such that T^ + e^ = 0. It follows from ( 18 1 that X^o + e 2 = 
0. By using (22), we can rewrite (23) as 



(c T A h iei + c T 4 +1 £ ) > a(c T A h 2 e 2 + c T A h 2 +^ ). 



(24) 



Define p L = cJ ' A\e\ + c t A^ +1 ^q for i = 1,2 and observe that there exists a > satisfying the above 
inequality unless (p\ ^ and p 2 > 0) or (pi < and p 2 ^ 0). However, neither of these two cases can occur 
due to the definition of £o and (16) with ^ = Co- 



5.1.2. Proof of Theorem {3J\B 



On the one hand, right-uniqueness of Filippov solutions (statement [lj necessitates that the implication 

±j = AiXi + ej, ccj(0) = £ £i(i) = ^(f) for alH ^ (25) 

holds for all (eWfnWj . On the other hand, it follows from statement [H] and Q that Wf = W$ = {£ | 
T 1 ' l+1 ^ + e^ +1 = 0}. As such, statement M immediately follows from fl25j ). 



Proof o/ Theorem \3im 

Note that both statements [5] and [6] are invariant under time-reversal. As such, it is enough to show that 
the forward Caratheodory property holds for every Filippov solution. 

It follows from statement [5] that W° := Wf = W 2 - We claim first that forward Caratheodory property 
holds for every Filippov solution with the initial state £ € W° . To this end, note that statement [6] implies 
that 

£ € W° + ei = A 2 £ + e 2 . (26) 



Note also that the implication 



?eW° =^ Ait + et e W° 



(27) 



readily holds from the very definition of the sets for i = 1,2. Together with ( 26 ), this invariance property 
yields 

±i = A^i + e u Xi(0) = £ =*> x x (t) = x 2 {t) for alUeR (28) 

for all £ G W°. Let x* satisfy x* = A\x* +e\ with a;*(0) = £. Also let x be a Filippov solution of the system 
with x(0) = £. Then, there exists a function A : R — > [0, 1] such that 



x(t) = X(t)[Aix(t) + e x ] + (1 - A(i)) [A 2 x(i) 



■ e 2 



for almost all i € 1. It follows from (28) that 



x* (t) = A(i) [^n* (t) + ej + (1 - A(t)) [A 2 x* (t) + e 2 ] 
for all f € R. Define A(t) := A(i)Ai + (1 - X(t))A 2 . Then, we get 



(29) 



(30) 



-( 



\x*(t) - x(t)\\ 2 ) = 2(x*(t) - x(t),x*(t) - x(t)) 

= 2(x*(t) - x(t),A{t)(x*(t) - x(t))) a\\x*(t) - x(t)\\ 2 
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where a := 2max{||AAi + (1 — A) ^2 1 | A G [0, 1]}. Since x*(0) — x(0) = 0, the last inequality readily implies 
that x(t) = x*{t) for all ( £ K. In other words, x* is the unique Filippov solution for the initial state £. It 
follows from ([27) that x* (t) G W° for all t G K. Hence, x* is a forward Caratheodory solution and statement 
[| holds for every Filippov solution with the initial state £ G W . 

Next, we claim that if x is a Filippov solution for the initial state £ and t* G K is such that x(t*) G W° 
then x(i) G W° for all ( e M. To see this, note that x(t) = x(t + t*) is a Filippov solution for the initial 
state x(t*) G W°. As such, the above argument yields x(t) — x(t — t*) G W° for all tel. 

Therefore, it remains to show that forward Caratheodory property holds for every Filippov solution x 
with the property that 

x(t) # W° (31) 

for all tgR. Let x be such a Filippov solution and let t* G K. If c T x(t*) + / 7^ 0, it follows from statement 
[H] and continuity of x that there exists e t » > such that Q holds. Suppose that c T x(t*) + f = 0. It follows 
from ( 31 ) and ( 12 ) that there exists an integer k with ^ q ^ h such that 

c T A[x{t* ) + c T A^ _1 ei = for alll = 1, 2, . . . , q (32) 
c T A« +1 a;(r)+c T ^ei t^O. (33) 

Now, suppose that 

c T A« +1 .T(r) + c T A? ei > 0. (34) 

It follows from statement [5] that 

c T A l 2 x(f) + c T A e 2 - 1 e 2 = for all I = 1, 2, . . . , q (35) 
c T A« +1 a;(r ) + c T A\e 2 > 0. (36) 

Since x is continuous, there must exist Cf > such that 

c T A\ +1 x{t) + c T A\e x > (37) 
c T A q 2 +1 x{t) + c T A\e 2 > (38) 

for all t G {t*,t* + e t .). 

For A G [0, 1], we define A(\) := XA X + (1 - X)A 2 , e(A) := Aei + (1 - A)e 2 . We also define 

Go = Ho := {/}, 
G k := {AiG' I G' G G k -i,i = 1,2} for k ^ 1, 
H fe := {A(A)iJ' I A G [0, 1] and i?' G H fe _i} for fc > 1. 

For each fc > 1 and (G, i? ) G Gfc x H fc of the form 

G = A ifc A,^ •••A il , i? = 4(A fc )A(A fc _i) • • • A(Xi), 

we define 

e G := c T A ik A lk _ 1 ■ ■ ■ A i2 e n and e H := c T A(X k )A(Xk-i) ■ ■ ■ A(A 2 )e(Ai) 

with the convention that e/ := / for k = 0. Finally, define Gk '■= {(G,eo) \ G G Gfc} and Hk '■= 
{(H,e H ) I H G H fc }. 

Note that Gfe C Hfc and C/fe C Hfe for all k ^ 0. Also note that conv(£/fc) C conv("Hfc) and that 
T-Lk C conv(C?fe). Hence, we have 

conv("H fe ) = conv(£fe) (39) 

for every k ^ 0. 
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Now, we claim that for each £ E {0, 1, . . . , q} and each GgGf 

c T Gx{t*) + e G = 0. (40) 

To prove this claim, we make an induction on I. The case £ = is evident. Suppose that the claim holds for 
all £ with ^ £ < p ^ q. Let G E G p . In the cases G = A\ and G = Af, the claim readily follows from (32 1 
and |55), respectively. Otherwise, G = A^A^G' for some i E {1, 2}, 2 sC fc x + fc 2 sC p, and G" € G p - kl - k2 . 
Then, it follows from statement [5] that there exists a positive number a such that 

c T Gx{t*) + e G = a(c T Gx{t*) + e G ) (41) 

where G — A^ 1+k2 G' . By repeating the same argument, one can prove the existence of a positive number 
a G such that 

c T Gx(t*) + e G = a(c T A p lX (t*) + jA^e,). (42) 



Then, the claim follows from either (|32j) or (|35j). The very same argument employed in the last step of the 
above induction yields 

c T Gx{t*) + e G > (43) 

for all G € Gq+i- Since x is continuous and G g +i is a finite set, we can conclude that there exists a positive 
number e* such that 

c T Gx{t) + e G > (44) 

for all G € Q q+ i and for all t e (t* , t* + e*). 
From (39), we further get 

c T Hx(t*) + e H = (45) 
for all H £ Mi and I = 0, 1, . . . , q, and also that 

c T Hx{t) +e H > (46) 

for all H e M q+1 and for all t £ (t* , t* + e*). 

c T Hx(t) E {c T H'x(t) + e H > | H' € H 9+1 }. (47) 



Let HeH,. Note that 



Also note that 



c T Hx{t) + e H = c T Hx(t* ) + e H + J c T Hx{s) ds. (48) 



Then, it follows from (|45j, (|46j, and (|47J that 

c T Hx{t) +e H > (49) 

for all t E (t* ,t* + e*). The very same argument can be repeated for H E Hf with £ = q — 1,...,1,0 to 
obtain 

c T Hx(t) + e H > (50) 
for all -ff E H^, £ = 0, 1, . . . , q, and t E (t*,t* + e*). In particular, we obtain 

c T x(i) + / > (51) 

for all f E (t*, i* + e*) with the choice of £ = 0. Therefore, x is a forward Caratheodory solution. 
The case 

c T A q 1 +1 x(t*) + c T A q 1 e 1 <0 (52) 
can be proven by using the above arguments in a similar fashion. 
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5.1.4- Proof of Theorem \3l\D 



As we have just shown, statements [5] and [6] imply that every Filippov solution is a forward Caratheodory 
solution. Then, it is enough to show that every forward Caratheodory solution is right-unique. 

Let x%,X2 be two forward Caratheodory solutions with a;i(0) = x 2 (0) = £. Then, there exists e > such 
that 

±i(t) = A lXl (t) + e h {-\) l [c T Xl {t) + /] > (53a) 
x 2 (t) = A 3 x 2 {t) + e h (-iy[c T x 2 (t) + f] > (53b) 

for all t £ [0, e). In case i = j, we have readily X\(t) = x 2 (t) for all t £ [0, e). In case i ^ j, we can assume, 
without loss generality, that i = 1 and j = 2. Then, it follows from (53) that £ £ Wi and £ £ W 2 - In other 
words, we have 



(-1)' 



for i £ {1,2}. Then, it follows from statement [5] that 



r c T - 




f 


c T Ai 




T 

c ei 




i + 




c T A h + 1 _ 




c T A>ie K 



y o 



(54) 



r c T - 




/ 


c T Ai 




T 

c ei 








P r A h i +1 _ 




_c T Ale K 



(55) 



for i £ {1, 2} and hence £ £ W° = W° = W 2 . As proven in 5.1.3 this means that xi(t) = x 2 {t) = x*(t) for 
all t £ R where x* satisfies x* — A x x* + e 1 with x*(0) = £. 

5.2. Proof of Theorem \3.4\ 

Since M is a lower triangular matrix with positive diagonal elements, it can be partitioned as 



M = 



Mi 

* m 33 



mix 








* 


m 22 





* 


* 





where ma > for i = 1, . . . , 3. 
It follows from ( 10 ) that either 



or 



c T A 1 _ 




e T A 2 _ 


) 


" / " 

T 

c e\ 


> Mi 


" / " 

T 

c e 2 


' c T ' 

c T A K 


= Mi 


' c T ' 

e T A 2 _ 


) 


' / " 

T 

c e\ 


= Mi 


' / " 

T 

c e 2 



(56) 



(57) 



holds. The claim follows readily for the former case due to Theorem |3.3| Therefore, it remains to prove the 
claim when (571 holds. 



Let £ £ R n \r2. Due to (10), (57) and the definition of f2, only the following cases are possible: 

1. c T £ + /^0. 

2. c T £ + / = 0, (e T A^ + c T ei)(c T v4 2 £ + c T e 2 ) > 0. 

3. c T £ + / = c T Ai£ + c T ei = c T A 2 £ + c T e 2 = and 

{c T A\Z + c T A iei )(c T AU + e T A 2 e 2 ) > 0. 
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4. c T £ + / = c T Ai£ + c T ei = c T A 2 £ + c T e 2 = and 

c 7 ^ + c T A iei > 0, + c T A 2 e 2 < 0. 



(58) 



By similar arguments to those in the proof of Theorem |3.1|C[ one can conclude the right-uniqueness of 
Filippov solutions for the initial state £ for the first three above-mentioned cases. 

For the last case, we claim that if a; is a Filippov solution with x(0) = £ then there exists e > such 
that 

c T x(t) + f = (59a) 
c T A lX (t) + c T ei = c T A 2 x{t) + c T e 2 = (59b) 
for all t £ [0, e). Note that the conditions (p39"| would imply that x satisfies the differential inclusion 



x £ < 



{Axx + ei} iic T x + f<0 

conv{Aix + ei,A 2 x + e 2 } if c T x + f = 
^{A 2 x + e 2 } if 5 T x + / > 



(60) 



on the interval [0, e) where c T = c T Ai and / = c T e\. As such, right- uniqueness of x would follow from (58 1 
and [1 Theorem 2.10.2] since c T £ + / = and c T A^ + c T ei > 0. 

Therefore, it is enough to prove the existence of e > such that the conditions (59) are satisfied. To do 
so, we first note that one can find positive numbers S and e such that 



c 1 A(x(t) + c 1 Aiei ^ 5, c 1 A 2 2 x{i) + c 1 A 2 e 2 sC -5 



(61) 



for all t £ [0, e) due to (58), x(0) = £ and the continuity of x. To show that the conditions (59) are satisfied 
with this choice of e > 0, consider the functions Vi,V 2 : K™ — > M defined by 



V 1 (z) = 



(c T z + f)(c T Alz + c T A iei ) ifc T z + /^0 
(c T z + f)(c T A 2 2 z + c T A 2 e 2 ) if c T z + f > 0, 



V 2 (z) = 



c T A 1 z + c T ei ifc T z + /<0 

m 22 (c T A 2 z + c T e 2 ) if c T z + f > 0. 



Also consider V(t) := -V|(ar(t)) - Fi(x(t)). From qsTJ, we get 

K(t) > 5|c T x(t) + /| + ^ 2 2 (x(<)) ^ 



(62) 



for all t £ [0,e). Note that the relation 

V(t) s$ AF(t) (63) 

for almost all t £ [0, e) and for some A > would imply V(t) = for all t £ [0, e) since V(0) = 0. Together 
with (62) and the definition of V 2 , this would imply that the conditions (59) are satisfied. As such, it is 
enough to show ( 63 ) in order to complete the proof. 

Define r := max{||a;(t)|| | t £ [0, e]} and B(0, r) :— {x | \\x\\ ^ r}. Note that V\ is readily continuous and 
also that V 2 is continuous due to ( [57] ). Furthermore, both V\ and V 2 are Lipschitz on B(0,r). Together with 
the absolute continuity of x, this implies that Vi(x(-)) and V 2 (x(-)) are even absolutely continuous on [0, e). 

Let A be the set of all t £ [0, e) for which x, Vi(x(-)), and V 2 (x(-)) are differentiable. Due to the absolute 
continuity of these functions, the set [0, e) \ A is of measure zero. Thus, it suffices to prove that (63) holds 
for all t £ A. 

To do so, let t* £ A. First, suppose that c T x(t*) + f ^ 0. In this case, we get 
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(64) 



where L := max{|c T Afx(t) + c T A 2 ei\ : t € [0, e],i = 1, 2}. In other words, ([63]) holds for t*. 
Now, suppose that c T x(t*) + f = 0. The Taylor expansion of x around t* yields 

c T x{t* + t) + f = Tc T x{t*) + c t o(t) 



(65) 



for all t sufficiently close to where lim T _ j . °( r )/ T = 0. It follows from (65) that the set of the one-sided 
derivatives of Vi (&(•)) at t* is equal to 



{c 1 x(t*)(c T A 2 x{t*) + c 1 A^) | * = 1,2}. 



(66) 



Since V\{x{-)) is differentiable at t*, this set must be a singleton. In view of (61), this can happen only if 
c T x(t*) = 0. This means that V\(x(t*)) — 0. Since x is a Filippov solution, we get 



€ conv{c T Axx(t*) + c T ei, c T A 2 x(t*) + c T e 2 }. 
By post-multiplying the first equation in ( p37"| by x(t*) and adding to the second, we obtain 

c T Aix(t*) + c T ei = m 2 2{c T A 2 x(t*) + c T e 2 ). 



(67) 



From (67) and (68), we get 

c T Aix(i*) + c T ei = c T A 2 x(t*) + c T e 2 = 0. 
From the definition of V 2 , this yields V 2 (x(t*)) = and further 

j t {V 2 2 {x{t*))) = 0, V{t*) = j t (Vi(x(t*))) - jVi{x{t*)) = 0. 



Clearly, (63) holds for t* . 



6. Conclusions 

We studied existence, uniqueness and nature of solutions Caratheodory and Filippov solutions for bimodal 
(possibly discontinuous) piecewise affine systems in a differential inclusion setting. First, we showed that the 
typical conditions that are employed in the context of general differential inclusions in order to guarantee 
uniqueness of Filippov are quite restrictive in the context of piecewise affine systems. Then, we presented a 
set of necessary and a set of sufficient conditions that ensure uniqueness of Filippov solutions for bimodal 
piecewise affine systems. By investigating the relationships between Caratheodory and Filippov under the 
presented condition, we provide conditions that rule our the so-called Zeno behavior. Possible extensions of 
the main results of this paper to general piecewise affine dynamical systems with external inputs emerge as 
future research directions. 
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